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Abstract 

We consider the spin-S ferromagnetic Heisenberg model in three dimensions, in the absence of an 
external field. Spin wave theory suggests that in a suitable temperature regime the system behaves 
effectively as a system of non-interacting bosons (magnons). We prove this fact at the level of the 
specific free energy: if S — > 00 and the inverse temperature /3 — > in such a way that j3S stays 
constant, we rigorously show that the free energy per unit volume converges to the one suggested by 
spin wave theory. The proof is based on the localization of the system in small boxes and on upper 
and lower bounds on the local free energy, and it also provides explicit error bounds on the remainder. 

1 Introduction 

An important open problem in theoretical and mathematical physics is the proof of long range order 
in the three-dimensional (3D) quantum Heisenberg ferromagnet (FM) at low temperatures. While 
the existence of long range order at low temperatures in the classical Heisenberg model and in the 
quantum Heisenberg antiferromagnet can be proved by reflection positivity methods [FSS, DLS], the 
broken phase of the quantum ferromagnet eluded any rigorous treatment so far. 

From a heuristic point of view, a very useful and suggestive representation of the quantum FM is in 
terms of spin waves, an idea first introduced by Bloch in his seminal work [Bl, B2]. The spin waves 
are the lowest energy excitations, which give the dominant contribution to the free energy at low 
temperatures; they satisfy a Bose statistics and are in many respects the analogues of the phonons in 
crystals (see, e.g., [Ke] for a classical and comprehensive review). Bloch's theory was later generalized 
in several directions by Herring and Kittel [HK], Holstein and Primakoff [HP], Dyson [Dl], and it was 
used, among other things, to compute a low temperature expansion for the spontaneous magnetization 
in zero external magnetic field: after a few erroneous attempts [Kr, Op, Sc, VK], Dyson's result [D2] 
was confirmed by a number of different methods [CH, KL, M, Og, RL, Sz, VLP, W, YW, Z] and 
further extended, more recently, by the effective Lagrangian method [Hi, H2]. The conclusion is that 
at low energies the corrections to the simple Bloch's theory coming from the interactions among spin 
waves are so small that for most practical purposes the linear theory is enough, both in the presence 
or in the absence of an external magnetic field. 

While physically Bloch's theory is accepted and in good agreement with experiments, from a more 
mathematical point of view there is no confirmation of its correctness yet. It is fair to say that the 
current mathematical methods of quantum many body systems are still far from allowing us to prove 
the existence of a spontaneous magnetization in the quantum Heisenberg ferromagnet and to possibly 
confirm the exactness of Dyson's computation [D2]. 

An easier, but still unsolved problem is to prove the correctness of spin wave theory at the level of 
thermodynamic rather than correlation functions, at least in the large spin limit, where the quantum 
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FM is (formally) close to its classical counterpart [L] or to the free Bose gas [HP, vHBW], depending 
on the temperature regime which we look the system at. Some progress in this direction is due to 
Conlon and Solovej [CS1], who proved that the free energy of the quantum Heisenberg FM in the 
presence of a large magnetic field is the same as the one of a free Bose gas, asymptotically for large on- 
site spin S. For small (i.e., 1/2) spin and small temperatures, Conlon and Solovej [CS2] and Toth [T] 
derived an upper bound on the free energy, close but not exactly equal to Bloch's theory prediction. 
In the large spin limit, it was unclear whether the rigorous justification of spin wave theory by Conlon 
and Solovej [CS1] could be extended to the case of zero magnetic field, which is much more subtle 
and unstable with respect to temperature-induced excitations, due to the global rotational symmetry 
of the model. For the classical counterpart of this problem, that is the convergence of the large-S 
classical Heisenberg model to the gaussian one, Conlon and Solovej managed to show the uniformity 
of their results as the magnetic field is taken to zero [CS3] (an even stronger convergence result to 
the gaussian model for the case of the classical rotator model is due to Bricmont et al [BFLLS]). 
However, their proof did not apply to the quantum case, see the comments after [CS3, Theorem 2.1]. 
The large spin limit for the XXZ chain in the absence of a magnetic spin was studied by Michoel and 
Nachtergaele [MN1, MN2], but their proof required the presence of an energy gap in the spectrum of 
the spin chain. 

In this paper, we solve the issue by extending these results to the 3D quantum rotationally 
invariant (gapless) case; namely, we show that the free energy of the 3D quantum Heisenberg FM at 
temperatures of order S converges as S — ¥ oo to the one of the non-interacting magnon gas. This 
rigorously proves (in a weak sense) the validity of the spin wave approximation at leading order in S 
for S large. We use the Holstein-Primakoff representation of the quantum Heisenberg FM in terms of 
interacting bosons, and we take advantage of the methods and the ideas recently developed by Lieb, 
Seiringer and Yngvason to treat the interacting Bose gas at low densities [LS, LSY, LY] (see [LSSY] 
for a review): key ingredients of our proof are: (i) a localization bound, allowing us to coarse grain 
the system in finite boxes; (ii) a priori bounds on the local energy, which make use of the energy gap 
in the box and allow us to drop the states with total spin far from the maximum. 

The paper is organized as follows. In Section 1.1 we define the model and state our main results. 
In Section 1.2 we review the Holstein-Primakoff representation and briefly review spin wave theory. 
In Section 2 we prove our main theorem by deriving suitable lower (Section 2.1) and upper (Section 
2.2) bounds. 

1.1 The Model and the Main Result 

We consider a ferromagnetic Heisenberg model with nearest neighbor interactions which is associated 
with the Hamiltonian 1 

Hr-.=j Yl (s 2 -Sx-s y ), (i.i) 

<x,y)cA 

where J > 0, A is a cubic portion of side length L of 1? with periodic boundary conditions, the sum 
runs over the nearest neighbor pairs in A, and S x is a three-components spin S operator. This means 
that the components of S x , denoted by 5 X , j = 1, 2, 3, satisfy: 

[Si, S$] = ie ]kl S l J^ y , S 2 X = (Si) 2 + (S 2 ) 2 + (S x ) 2 = S(S + 1) , (1.2) 

where in the first equation Ekim is the completely antisymmetric symbol, while in the second equation 
S is the size of the spin, with 25 an integer. We shall denote by M'k the Hilbert space of spin 
configurations in A such that S x = S X (S X + 1) with S x = 5"; a convenient basis for this space is 
({Sx}) := (g) xeA \S'i), with -S < S'i < S, Vx G A. Note that the Hamiltonian (1.1) is normalized 
in such a way that the ground state energy is zero. 

Our main object of interest is the free-energy per site 

f(S,P,A) :=-Jn log Z(S,P, A) , with £(S,/3,A) = Tr^ A (exp{-/?tfr r }), (1-3) 



1 Wc use the convention of distinguishing (when needed) between operators and numbers or eigenvalues by means of a 
on top of the first ones. As usual we also denote vectors by bold letters. 



Free Energy of the Hciscnbcrg Model - Correggi, Giuliani - July 18, 2012 



3 



and the thermodynamic limit of the specific free energy, namely, f(S,/3) := limA->oo f(S,/3,A). Our 
goal is to compute f(S,/3) to the leading order in 5* for 5* large in a suitable temperature regime (to 
be discussed below) and to prove that it coincides within explicitly estimated errors with the free 
energy of a free Bose gas. More precisely, our main result, which generalizes a similar result in [CS1, 
Theorem 4.2] proven for a non-zero (in fact huge, i.e., h = O(S)) magnetic field, can be stated as 
follows. 

Theorem 1.1 (Free energy asymptotics). 

Assume that j3 = fiS^ 1 as S — > oo for some (3 > constant independent of S. Then one has 

^--Jj^^ g (l-e-^) + o(s-^(lo g S)^), (1.4) 

B 

where e(k) := X/i=i (1 ~~ cos is the dispersion relation and B = [—ir, it] 3 is the first Brtllouin zone. 
Remark 1.1 (Zero temperature limit) 

By a direct inspection of the proof, it is easy to verify that all the estimates are uniform in j3 > 0, 
if /3 is well separated from 0. Hence the result can be straightforwardly extended to any ft — fi(S) 
function of S, provided that lim mfs-xx, /3(S) > 0, including the case that lims^oo = oo. This 
last case for the spin-1/2 Heisenberg ferromagnet was considered in [CS2] and [T], where only upper 
bounds on the free energy were proven. 

Remark 1.2 (External magnetic field) 

The addition of an external magnetic field would not affect the result proven in the Theorem above. 
In fact, as already mentioned, the case of a magnetic field h of order O(S) was already considered 
in [CS1]. However any magnetic field h <C S would produce a contribution to the r.h.s. of (1.4) of 
order h/S and thus would not affect the leading order term, irrespective of the sign of h. 

The intuition behind this result is based on a well-known bosonic representation of the Heisenberg 
model, first proposed by Holstein and Primakoff in [HP] , reviewed in the following section. 



1.2 Bose Gas Representation and Magnon Approximation 

It is well known since the pioneering work Holstain-Primakoff [HP] that the Heisenberg Hamiltonian 
can be rewritten in terms of suitable creation and annihilation operators, so that the excitations of 
the model can be described as Bose particles exactly as phonons in crystals. For any x € A we set 



1-: Hf' ^ =: V2S \I 1 ~ 4f iax ' S^=:aU^~S, (1.5) 

where aj., « x are bosonic creation and annihilation operators. We can associate with such Bose modes 
an Hilbert space isomorphic to the spin Hilbert space .Jfk by diagonalizing the number of particles 
operator n x := a^a x and we denote such a basis by |{n x }) := ® xeA Note that the constraint 

— S < < S translates into the requirement 

n x < 25*, (1.6) 

i.e., the Hilbert space is truncated and the occupation number per site can not exceed 2S. The 
isomorphism between the Hilbert spaces J%a and J?a can be implemented by means of the one-to-one 
correspondence \n x ) < — > 5 X = n x — S), so that, for instance, the ground state with all the spins 
pointing down is |{n x = 0}). 

The Hamiltonian H^ er can be rewritten in terms of the creation and annihilation operators as 

HT = JS £ \-ai(l-^) 1,2 (l-^) 1/ \ y + aU,-^aUWa y \, (1.7) 



x,yeA 

llx-yllA^ 1 
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where || ■ || is the euclidean distance on the torus A. For S large, this representation naturally leads 
to a decomposition of T-i p ^ r into "free" and "interacting" parts: 

nr=K: r A +icr, kj-.^sj e ^k-m, a.s) 

x,y6A 

ll x -y||A = i 

where rZ p ^ T is at least quartic in the Bose operators and is formally of relative size 1/S with respect 
to the quadratic part. The magnon approximation simply amounts to drop the interaction term K, p ^ r 
from the energy "H A er and study the free Bose gas so obtained. By means of the Fourier transform 
the free Hamiltonian 'Hq'a can De easily diagonalized: given any k £ A*, where A* denotes the dual 
lattice, i.e., k = ^-m, with < rrn < L, we set e(k) := 5^f =1 (l — cosfc) and define the creation 
and annihilation operators for the Fourier modes as fik := L~ 3//2 X^xeA ex P {*k ' x } fl x, so that the 
Hamiltonian (1.8) can be rewritten 'Hq'a = SJ X^keA* e (k)2] t ak. 

The free energy of such a Bose gas should be computed by taking the trace of exp{— PH pe ^} 
on the truncated Hilbert space J^a, where the occupation number of each site cannot exceed 25*. 
However, in the large S limit one expects that such a constraint can be removed in the calculation 
of the free-energy, up to higher order corrections, and in the thermodynamic limit this would yield 
exactly the dominant contribution in the r.h.s. of (1.4). 



2 Proofs 

In this section we prove Theorem 1.1. The asymptotics (1.4) is obtained by comparing suitable upper 
(Section 2.2) and lower (Section 2.1) bounds to the free energy of the ferromagnetic Heisenberg model. 



2.1 Lower Bound 

Here we prove the lower bound to the free energy. The main ingredients in the proof are a localization 
of the energy into Neumann boxes and two simple lower bounds on the localized Hamiltonian: the 
first one uses the energy gap in the box to show that the states with total spin far from the maximum 
have large energy; the second one is a rough bound on the contribution from the bosonic interaction 
part, which will be useful only after having restricted the trace to states with third component of the 
total spin close to the minimum value (the proofs of these two estimates are deferred to Appendix 
A). 

Proposition 2.1 (Free energy lower bound). 

Assume that P = PS^ 1 as S — > oo for some P > constant. Then one has 



f(S,[3) > 1 f dk 



-log (l - e -* J ««) - O (s- 1/8 (log5) 1/4 ) . (2.1) 



S ~ p J (2tt) 

B 

Proof. The first step towards the proof of (2.1) is a localization of the energy into small boxes of side 
length £ with Neumann conditions at the boundary. We partition the big box A into boxes A, of side 
I and use the positivity of the bond energy for all the bonds connecting a site in A, with a site in 
' 7^ ]• Correspondingly, we bound the original Hamiltonian H from below as H > where 
Hf[. depends only on the degrees of freedom associated with the spins in A* and has free conditions 
at the boundary: 

< = | E (S 2 -S x -S y ), (2.2) 

x,y£Ai 
x — y| — 1 

where | ■ | is the euclidean distance on Z 3 . The Hilbert space on which HZ is assume to act (i.e., the 
restriction of Jtfk to Ai) will be denoted by Jtfki- 

Obviously, all the Hamiltonians Hj^. commute among each other. Therefore, 

Z(S,P,A)< [Tr^r Ai {e^ V [-pHlY)] L3/l3 ^.Z N {S,P,k 1 ) L3 / e3 , (2.3) 
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where the trace Tr^r Ai is only over the spin degrees of freedom within Ai. 

In the computation of Trj^ Ai exp{— fiH^ } we now distinguish between the states with total spin 
close to the maximum, from those with "small spin", which have a big energy and, thus, give a small 
contribution to the free energy. Note that the Hamiltonian is invariant under global rotations, 
that is commutes with the three components of the total spin St := X^xgAi ^ x m Ai. Therefore, 
is block diagonal with respect to the decomposition J^. 1 = (Bs^=o©s 3 -_s J^s T sf.> where 
J^s T sf, i s ^ ne subspace of Jij on which S T = St{St + 1) and S T = S T . On each M' St S 3 , the 
Hamiltonian can be bounded from below as follows, independently of the value of S T (see Appendix 
A for a proof). 

Proposition 2.2 (Lower bound on H^^. 

There exists a positive constant c > such that, if S — £~ s St > 3Jc~ x £ 2 , then 

<L SrS3 ^ rf5 ( 5 -fr)- ( 2 - 4 ) 

We now split the trace of interest in two parts: 

Z N {p,SM)= E E Tr ^ Sr , S 3 exp{-0flft}+tt, (2.5) 

S T >S*£ 3 S 3 =-S T 

where S 1 * < S is a parameter which is going to be fixed later, and the rest TZ can be bounded by 
means of (2.4) as follows: 



a 

TZ = 

s T <s„« 3 s 3 =-s. 



E E Tr^ STs| exp{-/3<}<(25+lfe-^ S ( 1 -^), (2.6) 



provided that S — S, > 3Jc 1 £ 2 , which can be satisfied by picking 

5, _ 2l 2 log(2S+l) 

S - CPS ( ] 

and £, S large enough, so that TZ < (2S + 1)"^ « 0(1). 

We are now left with the main contribution to the trace (the one involving St > S f £ 3 ), to 
be called Z(j3, S, Ai). By using once again the fact that commutes with St, we find that 

Tr jz> s s3 exp { — /3i?Ai } i s independent 2 of S T , so that 

si 3 

Z09,S,Ai)= E ( 2 ^ + l)^ Sr ,_ ST exp (2.8) 
s T =s^ 3 

We can now apply the boson representation given in Eq.(1.5), which implies 

Ztfi, S,Ai)< {2S£ 3 + 1) E (i n *l | ex P {-Wm } | {"»}) , (2.9) 

{?i x },xgA 1 
£n x <(S-S*)f 3 



2 This can be proved as follows. Let us indicate by xi, . . . , x f 3 the sites of Ai labeled in lexicographic order. By the 
theory of the composition of angular momenta, a bona fide basis for Jf?^ is provided by the common eigenvectors of 
(S X1 + S X2 ) 2 , (S X1 + S X2 + S X3 ) 2 , . . ., (S X1 + • ■ ■ + S Xf3 i ) 2 , S 2 , Sf. In other words, the eigenvalues of (S X1 + S X2 ) 2 , 

(S X1 + S X2 + S X3 ) 2 , . . ., (Sxi + • • • + Sx,3 j) 2 can be used as good quantum numbers for classifying the states of S 3 . 
Note that the operators associated with these quantum numbers are all scalars, i.e., they commute with the three components 
of St: therefore, the eigenvectors of on J%g T g3 are invariant under the action of St, which implies in particular that 

Tr jg> s s3 exp{— /3iZ^} is independent of S'^.. 
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where Ha 1 is the bosonic hamiltonian in Ai with Neumann boundary conditions (i.e., it is given in 
Eq. (1.7) with A replaced by Ai and the condition ||x — y || a = 1 replaced by |x — y| = 1). Inspired 
byEqs. (1.7)-(1.8), we rewrite H\ t = ?< Al + /C£ , with H^ Al = JS £ <x , y)cAl ( a x - 4)(«* - o y ) 
and X3 Al the interaction part, which can be bounded as follows (see Appendix A for a proof). 

Proposition 2.3 (Estimate of /C Al ). 
There exists a finite constant C such that 



\K.AA<CNi^ 



(2.10) 



where AT Al ■= X] 



x6Ai 



Using this estimate in (2.9) together with the fact that N Al < (8- S,)l 3 = (2/c/3)£ 5 log(2S + 1), 
we get for a suitable constant C 



Z(/3,S, Ai) < (2ST 5 + l)e 



C'f 10 S _1 log 2 s 



Zo(J3,S,Ai), 



(2.11) 



where Zo(/3, S, Ai) stands for the partition functions of a free Bose gas in a box Ai with free conditions 
at the boundary and constraint on the total number of particles AT Al < (S — S*) £ 3 . 

The estimate of Zq is very simple: we rewrite the partition function in terms of the Neumann 
Fourier modes k G A* N (see Appendix A) and drop the constraint on the total particle number 
except for the mode k = and obtain 



Z (p,S,Ai)< Y, exp|-/3J Y eOc)h k } = (S-S ir )t 
h <{S-s„)e 3 

Putting together all the estimates, we obtain 



n 

keAJ N 



l-exp{-/3e(k)} 

(2.12) 



/(^,A )> ^ E 



(const. 



f^ogS) 2 + log(£ 3 S) 



(2.13) 



Finally, replacing the Riemann sum in the r.h.s. by the corresponding integral, we get 



ms,A) > i 



dk 



3^(1 



-y3Js(k) 



const. 



f(\o g S) 2 \og(£ 3 S) 1 
S P I 



(2.14) 



Optimizing over I yields I = S ,1//8 (log S) and the desired lower bound on the specific free energy 



□ 



2.2 Upper Bound 

In this section we complete the proof of Theorem 1.1 by deriving a suitable (rough) upper bound 
to the free energy. The most relevant steps in the proof are an energy localization into Dirichlet 
boxes, the introduction of a fictitious magnetic field (which will be removed at the end) to take into 
account the constraint on the bosonic Hilbert space and the estimate proven in Proposition 2.3 on 
the interaction. A finer upper bound, supposedly optimal up to corrections of the order 1/S included, 
will be presented elsewhere [CGS]. 

Proposition 2.4 (Free energy upper bound). 

Assume that /3 = as S — > oo for some (3 > constant. Then one has 

^ <~f ^lo e (l-e-^)+o(s-^(lo g Sf/ 3 ) . (2.15) 

8 
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Proof. We start by adding an external magnetic field h > to the Hamiltonian: exploiting the 
positivity of S-J + S, one has the trivial inequality 

Z(fi,S,A) > Z h {t3,S,k) := Tr,*» A (exp {-/?#,,}) (2.16) 

for any h > 0, where H h = H + hJ2 xeA (^i + We fix £ G N and define the corridor C of width 
1 as the minimal connected set on A that contains {x G A : x = n£, n 6 Z 3 } (here we assume for 
simplicity that L is divisible by £). Note that A \ C is a union of boxes Ai of side length £ — 1. We 
localize the energy in the boxes Ai with Dirichlet boundary conditions, by proceeding as follows. Let 
us denote by -Xc and J#a\c the Hilbert spaces generated by the spins in C and A \ C, respectively. 
Note that, if |4-c) is the state of ,_?t?c such that = — S, Vx G C, and l^} is a generic state of ^"a\c, 
then ((Id ® W)H h (\lc) ® |^» = <V>I E< #£a, IV>>, where 



#£ Ai = + /» + S ) + J E + S ^ , (2.17) 

xgA, xGSA; 

and the first sum runs over pairs of nearest neighbor sites, both belonging to Ai, while the second 
term is the boundary contribution, which should be thought of as an interaction term between the 
sites at the boundary of Ai with the neighboring sites in the corridor C. Obviously, the Hamiltonians 
H®. all commute among each other and can be naturally thought as operators on j#a ; . We also 

denote by Ef and \Ef) the eigenvalues and eigenvectors of H®. and set |{_Ef}) := (£)f = ( e 
Given these definitions, we note that 

Z h (P,S,A) > J2 ((U®({E?}\)e- 0H »(\lc)®\{EP})) 
{eP} 

> e*p{-KM®({E?}\)M\k)®\{E?}))} 

{BP} 



I? I 

E 



where, in order to go from the first to the second line, we used Jensen inequality to lift the expectation 
value to the exponent, while to go from the second to the third we used that the expression in braces 
is equal to -j3 ({E?}\ fl£ \{E?}). 

Now we use the bosonic representation, drop from the partition function Z®(f3, S, Ai) the con- 
tribution from large total occupation number A r A 1 , estimate the interaction IC Al and finally re- 
store the missing part of the partition function. To this purpose we introduce a new parame- 
ter N <C S, which is going to be chosen later and notice that the bosonic analogue of H h Al is 
T^Ai = Wa 1 + ^ExeA! n x + •S'^X]xeaA 1 nx - Now, in analogy with the case of Neumann boundary 
conditions, we can rewrite = W^h Ai + K*hAi > wnere Hoh Ai stands for the quadratic part 

of Hhhx > while /Cj^Ai is a t least quartic in the Bose operators and can be bounded exactly as in 
Eq.(2.10). Therefore, 



2^(^S,Ax)>e X p^-^^^Tr, 



.( Y, ™x< N) exp{-^^, Al } 



(2.18) 



where is the unconstrained bosonic Hilbert space on Ai, i.e. the Bose particle number n x is 
unbounded on : note that the removal of the constraint is irrelevant, because N <C S and, there- 
fore, the condition that n x < 25* is automatically satisfied. Now we rewrite 1 (£ x£Al n * — N) = 
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1 — 1 (^xeAi nx > Nj and bound the contribution from the states with more than N particles as 
y~] ({n*} expj-jSH^Ajj I {n*}) < ^ ({>x}j exp { - fiS^h ^ n x |j{n x } 



x6Ai 



< e 2 



PS~ L hN 



n e f 

xgA! n x >0 



-PS~ L hn, 



<e-&*~ 



(2.19) 



provided that ft/S < 1 and S~ 1 hN > f 3 log(Sy/i). The estimate (2.18) thus becomes 
^(/3,S J A 1 )> e - C ^ 2/S (l-e-^ S "^)Tr jr o i (exp{-^, Al }) . 



(2.20) 



In order to compute Trjr Ai (exp {— flHo,h,Ai }) we S° to Fourier space: assuming for defmiteness 
that Ai = {x £ 1? : Xi = 1, . . . ,i — 1}, we rewrite H o h ^ by applying the Fourier transform 
a k = ExeAj V 3 k(x)a x : , where y>k(x) = H i=1 <pk t (as»), ^k(aO = [2/(^— 1)] 1/2 sin(fc:r) and the set of 
momenta A* iZ3 is A 1£) = {fn, ngZ 3 : n» = 1, . . . ,£ — l}. After Fourier transform we get 



keAT 



r ttk 



(2.21) 



so that 



k)+h/SJ) 



(2.22) 



Putting all the estimates together gives 
logi^CS, S,Ai) < 

| Ai| 



-0J(e(k)H 



">) + (const.) |L - i log (l - e-^ 5 " 1 ' 



kGA*,k^0 

dk 



(27T) 



^ log ( 1 — e 



-/SJe(k) 



+ (const.) 



PS P g V 



h 1 

5 + £ 



<(2.23) 



(2.24) 



To satisfy the condition 5 1 hN 3> £ 3 log(S/h) met during the proof we can pick h = N' 1 S£ 3 (log 5") 2 
Optimizing over A' and £ yields 

N = CS 2/3 (logs')" 2/3 , £ = CS 1/6 (logS)" 2/3 



/i = CS 5/6 (logSf"«S, (2.25) 



n2/3 



for a suitable constant C. The corresponding relative error term in the energy is 0(S log 2 / 3 S), 
which proves the proposition. □ 



A Two Technical Bounds 

In this appendix we prove Propositions 2.2 and 2.3. 

Proof of Proposition 2.2. Let us assume for defmiteness that Ai = {x 6 Z 3 : Xi = 1, . . . ,£}. We 
rewrite by applying the Fourier transform <S k = I] xgAl 0k(x)S x , where 0k(x) = nLi ^{xi), 

{l" 1 / 2 if it = 

(2/£)^ 2 cos(fc(x-I)), iik^Q, {AA) 

and the set of momenta A lN is 

Al N = {£n,ne Z 3 : = 0, ...,£- l} . (A.2) 
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Then a straightforward computation shows that, if e(k) = X^i=i( i — cosfo), 



= -3SJ(£ 3 — I 2 ) + J £ ( k )«5k • S k . (A.3) 



Now, for every k/0, one has e(k) > cq£ 2 , for a suitable Co > 0. Therefore, 

r w| ^ oot«3 i t c o c ocr*3 , r C0 



#aJ > -35J^ 3 + J-| V «S k • 5 k = -3SJ£ 3 + J- 2 

e 2 



S(S + l)t-^(S T + l) 



, (A.4) 



where in the last inequality we used PlancherePs identity J^keA* = X^xeAi ^x — ^(S + l)^ 3 

and the definition of total spin operator St = £ 3 ^ 2 So, from which 5o = £~ 3 St(St + 1) on 3tfs T - 
Now, the expression in square brackets in the r.h.s. of Eq. (A.4) can be bounded from below by 
S£ 3 (S - S T /£ 3 ), so that 

flf, I > -3SJ£ 3 + c JSe ( S - ^ j , (A.5) 

Therefore, for all St such that S - ^ 3 St > 6c,7^ 2 , we get ff£ I > %J£S(S - ^~ 3 5 T ), which 
proves the proposition with c := \Jcq. □ 
Proof of Proposition 2.3. We start by investigating the part of associated with the square roots 

in the first term in (1.7). Using the fact that A x , y := 1 — \J 1 — §§^\A — is a non-negative operator 
on the bosonic Hilbert space of interest, we get: 

2Sj\ J2 a x-4x,ya y | = 2Sj] a i A l/,y A l/,y a y \ < 2 $J a x-4x, y a x < 

(x,y)CA 1 (x.y)CAi (x,y) CAi 



< J Y, a x("-x + n y )a x < (const.) ( re x ) , (A. 6) 

(x.y)CAi ^xeA x ' 

where in the first line we used the Cauchy-Schwarz inequality, while to go from the first to the second 
line we used the simple fact that 1 — \/\ — x < x, Vx G [0, 1]. Finally, the remaining term in 

can be bounded trivially as |X)<x,y)cAi n x"-y| < (const.) ^J^xeAi n^j , which implies the desired 
estimate. □ 
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